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Preface 


The  study  reported  herein  was  conducted  by  the  U.S.  Amiy  Engineer  Waterways  Experiment  Station 
(WES)  during  fiscal  years  1991-1992  as  pan  of  the  Environmental  Characterization  for  Target 
Acquisition  Program.  This  program  was  jointly  funded  by  the  U.S.  Army  Aviation  Applied 
Technology  Directorate  (AATD),  Fort  Eustis,  VA,  and  by  Headquarters.  U.S.  Army  Corps  of 
Engineers,  Washington,  DC.  Mr.  Nyle  Wilcocks  was  the  AATD  Technical  Monitor. 

These  efforts  were  under  the  general  supervision  of  Dr.  John  Harrison,  Director,  Environmental 
Laboratory  (EL),  Dr.  Victor  Barber,  Acting  Chief,  Environmental  Systems  Division  (ESD),  EL,  and 
Mr.  H.  Wade  West,  Chief,  Environmental  Analysis  Group  (EAG),  ESD.  Mr.  Bruce  Sabol,  EAG,  was 
the  Principal  Investigator  for  this  project.  Dr.  John  O.  Curtis  and  Mr.  Sabol.  EAG,  prepared  this 
report,  liie  labor-intensive  task  of  assembling  Appendix  C  was  capably  performed  by  Mr.  Sean 
Brewer  (EAG  contraa  student). 

At  the  time  of  publication  of  this  report,  EHrector  of  WES  was  Dr.  Robert  W.  Whalin.  Commander 
was  COL  Leonard  G.  Hassell,  EN. 

This  report  should  be  cited  as  foUows: 

Curtis,  J.  0-,  and  Sabol,  B.  M.  (1993),  "Enviiotunental  characterization  for  target  acquisition; 

Rqx)tt  3,  New  concepts  for  evaluating  low-grazing  angle  radar  measurements,"  Technical 

Report  EL-93-9,  U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


Conversion  Factors,  Non-SI  to  SI  Units 
of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  of  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

meters 
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1  Introduction 


Background 


One  possible  scenario  for  modem  military  helicopter  operations  involves  the  use  of  both  natural  and 
man-made  ground  features  to  provide  cover  for  the  aircraft  until  the  pilot  decides  to  pop  up  and 
attempts  to  locate  and  fire  on  a  suspected  target.  Target  acquisition  and  weapon  delivery  procedures 
will  ultimately  be  automated,  making  use  of  a  number  of  different  sensors  whose  data  may  be  fused  to 
provide  the  greatest  opportimity  for  success.  The  design  of  such  sensors  and  sensor  fusion  algorithms 
can  be  optimized  by  a  thorough  understanding  of  terrain  and  target  interactions  with  each  other  and 
with  each  sensor  under  various  environmental  conditions. 

To  this  end,  the  U.S.  Anny  Aviation  Applied  Technology  Directorate  (AATD),  located  at 
Fort  Eustis,  VA.  managed  a  field  demonstration  of  multiple  sensors  called  the  Multi-Sensor  Fusion 
Demonstration  (MSFD)  whose  purpose  was  to  collea  target  acquisition  perfcimance  data  on  targets  on 
the  ground  using  data  from  individual  sensors  as  well  as  combinations  of  sensors.  The  MSFD  took 
place  in  February  and  March  1988  at  Fort  Hunter  Liggett,  CA.  Sensors  were  provided  and  tested  by 
two  oiganizations,  Martin  Marietta  Oriando  Aerospace  and  a  combined  team  from  Hughes  Aircraft 
Company  and  Texas  Instruments.  The  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
was  initially  tasked  by  AATD  to  provide  ground  truth  information  for  these  tests  and  to  conduct 
studies  on  the  impaa  of  environmental  conditions  on  sensor  performance,  particularly  within  the  visual 
and  thermal  infn^  portions  of  the  spectrum. 


Objectives  of  This  Study 

After  completion  of  the  MSFD,  WES  was  additionafty  tasked  with  providing  input  to  AATD  on 
how  test  sites  might  be  characterized  with  respea  to  sensors  that  operate  in  the  millimeter  wave  por¬ 
tion  of  the  spectrum.  The  WES  staff  felt  that  the  first  step  that  should  be  taken  to  respond  to  this  last 
uteking  was  to  closely  reexamine  existing  data  collected  by  millimeter  wave  sensors  with  the  goal  of 
assessing  the  impact  of  terrain  conditions  on  the  sensors’  ability  to  locate  and  identify  targets  of  inter¬ 
est  and  not  simply  measuring  their  aWlity  to  find  those  targets.  In  other  words,  the  objectives  of  this 
study  are  to  determine  whether  or  not  the  terrain  presents  target-like  features  to  the  sensor  and  to 
determine  if  tliat  behavior  can  be  quantified  and  modeled  so  that  one  might  be  able  to  predict  future 
test  and/or  battlefield  performance  for  a  given  millimeter  wave  sensor. 
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WES  has  devcioped  an  image  metrics  approach  to  characicn/ing  test  sites  within  the  visible  and 
thermal  infrared  portions  of  tlit  spectrum.'  Because  radar  data  like  that  collected  during  the  MSFD 
can  be  displayed  in  an  "image’'  form  (slant  range  versus  azimuUi  position),  the  most  logical  approach 
to  quantifying  laiget-like  features  within  these  images  was  to  apply  the  appropriate  target  feature 
metrics  (filters)  to  the  radar  data  collected  at  Fort  Hunter  Liggett.  Therefore,  one  path  taken  in  this 
study  was  to  acquire  as  much  MSFD  r.car  data  as  possible,  to  convert  those  data  to  "images,"  and  to 
apply  the  simplest  radar  target  featu'  metrics  to  those  images  to  quantify  target-like  features  within 
the  enviromnent. 

The  second  path  taken  to  characterize  test  site  terrain  with  respect  to  radar  scasors  was  to  use  a 
very  simple  reflectance  model  to  predict  the  back,scatier  response  of  radar  energy  that  is  due  only  to 
terrain  surfa'  v  geometry  at  the  test  site.  A  point  light  source  was  placed  at  the  two  locations  for  the 
MSFD  sensor  test-beds  and  a  simple  Lambertian  scattering  model  was  used  to  calculate  backscattered 
power  from  a  three-dimensional  finite  element  representation  of  the  test  site  terrain  surface.  These 
predictions  were  compared  qualitatively  with  real  data  to  determine  the  relative  impact  of  terrain 
geometry  on  radar  returns.  Furthermore,  where  available,  overlays  of  the  locatioas  of  trees  and  signifi¬ 
cant  shrubbery  were  prepared  to  supplement  the  qualitative  assessment  of  the  impact  of  environmental 
conditions  on  radar  data. 


Scope  of  Report 

Chapter  2  contains  a  description  of  the  radar  data  from  the  MSFD  that  were  acquired  by  WES  for 
this  study  as  well  as  some  details  about  how  those  data  had  to  be  processed  to  produce  "images."  A 
description  of  the  target  feature  metrics  used  in  this  snidy  are  found  in  (Thapter  3  along  with  the  results 
of  applying  those  filters  to  one  of  the  MSFD  tests.  Chapter  4  discusses  the  light  model  and  vegetation 
overlay  studies,  and  Ch^ter  5  summarizes  the  results  and  makes  recommendations  for  fiiture  studies 
and  data  collection  methodology. 

Three  appendices  are  included.  Appendix  A  contains  a  summary  of  statistics  associated  with 
applying  the  target  feature  metrics  to  all  of  the  available  radar  scenes.  Appendix  B  summarizes  target 
feature  metric  statistics  as  they  relate  to  the  man-made  targets  within  each  scene.  Appendix  C  is  a 
graphical  output  summary  that  contains  scene  filter  results,  terrain  contour  maps,  backscaner  prediction 
results,  and  vegetation  overlays  for  all  of  the  scenes. 


'  B.  M.  Sabol  and  S.  Rivera.  (1993).  "Environmental  characierizaiion  for  target  acquisition;  Report  2.  Analysis  of  EO 
imagery  (in  preparation),"  U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg.  MS. 
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2  MSFD  Radar  "images" 


Data  Supplied  to  WES 

Through  the  efforts  of  the  AATD,  radar  data  for  31  separate  MSFD  tests  were  supplied  to  WES  by 
Martin  Marietta.  Of  these  31  files,  which  represented  all  of  the  MSFD  data  possessed  by  Martin 
Marietta,  only  24  could  be  read  on  WES  data  processing  systems  without  errors.  Data  were  not  avail¬ 
able  from  the  Hughes/Texas  Instrument  team.  Details  on  the  format  of  the  Martin  data  and  the  proce¬ 
dures  required  to  generate  radar  "images"  from  these  data  are  presented  in  a  later  section. 

When  reduced  to  image  format,  the  Martin  radar  can  be  thought  of  as  illuminating  a  patch  of  terrain 
roughly  420  m  in  range  depth  and  encompassing  an  azimuth  sweep  of  about  24  deg.'  Table  1  con¬ 
tains  a  listing  of  all  of  the  useful  Martin  data  by  MSFD  test  numbers.  Also  contained  in  that  table  is 
an  average  range  and  an  azimuth  angle  (clockwise  from  north)  that  can  be  used  to  approximately 
locate  the  center  of  the  illuminated  area.  The  "High"  and  "Low"  site  designations  on  this  table  signify 
the  location  of  the  radar  system  for  each  test.  Depression  angles  from  the  radar  locations  to  the  center 
of  the  test  sites  never  exceeded  6  deg. 

Figures  1  and  2  contain  ground  plane  representations  of  the  approximate  area  covered  by  each  of 
the  MSFD  tests  for  the  radar  located  at  the  high  site  and  low  site,  respectively.  The  elevation  contour 
lines  superimposed  on  dtese  figures  are  at  2(X)-fi  intervals. 


Creation  of  Radar  "Images” 

The  31  MSFD  test  files  sent  to  WES  by  Martin  Marietta  consisted  of  range-corrected  (but  not  cali¬ 
brated)  injrftase  and  quadrature  frequency  domain  data  for  the  K,-band  circularly  polarized  radar 
system  us«l  at  the  MSFD.  Data  files  ranged  in  si2»  ftom  23  MBytes  to  73  MB)hes  and  contained 
multiple  azimuth  sweeps  of  each  test  area  and  usually  two  sets  of  different  range  gate  data.  The  term 
"range  gate"  as  used  in  this  report  refers  to  a  set  of  data  that  is  collected  by  the  receiver  and  recorder 
within  a  very  specific  range  of  time  delays  following  tlw  transmission  of  the  wave  form.  Each  range 
gate  contained  20  coarse  range  cells  that  were  defined  by  the  pulse  width  of  the  radar  wave  form. 
Furthermore,  the  data  were  collected  by  stepping  frequencies  in  both  an  up  and  down  direction,  the 
total  bandwidth  of  the  sweep  determining  a  theoretical  limit  on  range  resolution.  Data  were  also  col¬ 
lected  with  the  center  of  the  beam  at  two  different  elevation  settings. 


'  A  uble  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (metric)  units  is  presented  on  page  vi. 
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Table  1 

MSFO  Test  Data  Used  for  These  Studies 


Test  No 

Radar  Site 

Avg  Range 
m 

Central  Az 

Angle,  deg 

Ceniral  El 

Angle,  deg 

MS5502 

HIGH 

4120 

161  6 

2  83 

MS5503 

HIGH 

1780 

232  5 

5  68 

MS5504 

HIGH 

2460 

192  2 

4  78 

MS5601 

HIGH 

3015 

238  1 

3  63 

MS5602 

HIGH 

3960 

i;i  8 

2  97 

MS5603 

HIGH 

3650 

2429 

3,17 

MS5604 

HIGH 

3040 

174  7 

3  72 

MS5702 

HIGH 

4460 

161  9 

2  70 

MSS  703 

HIGH 

3936 

251  7 

2  77 

MS5704 

HIGH 

2340 

202,4 

4  75 

Msseoi 

HIGH 

2430 

2327 

3  58 

MS5802 

HIGH 

4460 

162  1 

2.72 

MS5603 

HIGH 

1450 

262.3 

5  50 

MS5804 

HIGH 

4450 

162.1 

2.20 

MS6001 

LOW 

3880 

284.7 

023 

MS6002 

LOW 

3995 

286  4 

0.33 

MS6004 

LOW 

3840 

285.8 

0.32 

MS6101 

LOW 

3960 

285.6 

027 

MS6102 

LOW 

3920 

286.5 

0.27 

MS6301S' 

LOW 

1960 

277.1 

0  75 

MS6301L' 

LOW 

3050 

277.1 

0.42 

MS6403 

LOW 

1555 

297.1 

0  70 

MS6603 

LOW 

2540 

2579 

0.37 

MS67C1 

LOW 

2440 

275.8 

067 

’  S  =  short  range;  L  = 

long  range. 
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Figure  1 .  MSFD  radar  test  high  site  data  c»verage 
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The  point  of  listing  all  of  these  variables  that  were  a  part  of  the  data  coUeaion  process  is  to 
emphasize  that  each  data  file  contained  an  enonnous  amount  of  useful  information.  For  the  purposes 
of  this  study,  much  of  those  data  were  redundant.  This  is  not  to  say  that  sophisticated  data  processing 
algorithms  could  not  make  good  use  of  multi^e  terrain  sweeps,  overlapping  elevations,  and  over¬ 
lapping  range  gates,  but  the  study  repotted  herein  was  not  intended  to  be  a  target  detection  optimiza¬ 
tion  effort.  Before  any  of  the  target  feature  metrics  that  form  the  heart  of  this  site  characterization 
procedure  could  be  applied  to  a  given  test,  some  kind  of  data  compression  had  to  be  perfonned. 
Acquiring  a  clear  understanding  of  how  the  data  were  arranged  and  developing  a  procedure  for  pro¬ 
ducing  a  smaller  useful  data  set  was  a  nontrivial  matter. 

Martin  Marietta  provided  WES  with  FORTRAN  codes  that  could  be  used  immediately  to  read 
record  headers  for  these  data  files  but  that  lacked  the  proprietary  subroutine  that  performed  the  inverse 
Fourier  Transforms  for  converting  to  the  time  (or  range)  domain.  In  addition,  it  was  not  discovered 
until  much  later  that  the  way  in  which  the  data  are  stored  in  the  frequency  domain  resulted  in  a  wrap¬ 
around  of  data  in  range;  i.e.,  one  had  to  know  the  algorithm  for  locating  the  center  of  each  coarse 
range  cell.  This  algorithm  was  finally  obtained  from  a  classified  Martin  Marietta  report. 

Without  going  into  any  great  detail,  the  basic  procedure  for  compressing  the  Martin  files  into  some¬ 
thing  more  useful  for  these  studies  was  to  create  a  fine-range  resolution,  two-dimensional  map  for  each 
radar  elevation  setting  and  for  each  test.  The  fine-range  map  was  established  by  first  identifying  a  pair 
of  azimuth  swee{»,  one  in  each  of  two  successive  range  gates,  with  azimuth  marks  on  the  record 
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headers  ihai  were  most  compatible.  This  was  done  by  reading  all  of  the  headers  and  computing  the 
minimum  sum  of  squares  of  differences  between  successive  range  gate  sweep  azimuth  marks;  in  other 
words,  find  the  pair  of  azimuth  sweeps  in  successive  range  gates  that  best  line  up.  An  inverse  Fast 
Fourier  Transform  routine  was  written  to  produce  fine-range  resolution  data  from  each  coarse  range 
cell,  overlapping  data  were  eliminated,  and  a  continuous  string  of  fine-range  results  was  produced  for 
each  azimuth  setting.  The  end  result  for  each  radar  elevation  or  depression  angle  setting  (the  elevation 
designated  "barO"  is  a  depression  angle  greater  than  the  elevation  designated  as  "barl")  was  a  multi¬ 
dimensional  airay  of  data  measurements  that  are  proportional  to  power  returned  to  the  radar  from  the 
terrain  and/or  targets.  There  were  1,400  fine-range  cells  within  each  of  33  different  azimuth  settings, 
with  a  separate  "image”  for  each  of  the  two  radar  polarization  combinations  (right  circular  received 
and  left  circular  received  for  right  circular  transmitted  signals).  In  this  way,  data  files  that  consisted  of 
tens  of  MBytes  of  data  were  compressed  to  two  files  (two  elevations)  of  about  370  kBytes  each  that 
could  readily  be  turned  into  two-dimensional  "images"  in  range-azimuth  space. 

The  intent  of  these  studies  is  not  to  find  a  way  to  improve  signal-to-noise  ratios  to  beaer  detect 
targets  within  natural  backgrounds,  but  to  answer  the  following  question;  Given  a  single  sweep  of  the 
terrain,  how  many  target-like  features  are  out  there?  Therefore,  an  attempt  was  not  made  to  optimize 
the  data  set  beyond  selecting  the  best  azimuth  alignment  of  successive  range  gates.  Although  not 
tested  during  this  exercise,  a  simple  averaging  of  the  multiple  sweeps  (subject  to  some  kind  of  azimuth 
alignment  reasoning)  is  highly  likely  to  greatly  improve  the  quality  of  data  for  stationary  targets. 

Some  kind  of  comparison  of  data  from  the  two  radar  elevation  settings  probably  could  also  be  used  to 
minimize  terrain  background  response  and  to  enhance  man-made  target  response. 

Numbers  alone  do  not  have  the  impact  of  graphics  in  transmitting  information  to  the  user.  There¬ 
fore,  several  FORTRAN  routines  were  written  not  only  to  j^ply  metrics  to  these  compressed  data  sets 
but  also  to  display  data  and  results  of  analyses  in  both  two-  and  three-dimensional  formats.  All  calcu¬ 
lations  for  these  studies  were  performed  on  a  Micro  VAX  II  630QE  computer  system,  while  color 
graphics  were  generated  on  a  Raster  Technologies  Model  One/360  graphics  display  system. 

Naturally,  in  range-azimuth  space,  each  fine-range  cell  will  be  represented  by  a  circular  arc,  and  the 
entire  test  data  set  will  be  represented  by  circular  segments  like  those  shown  on  Figure  1.  However, 
because  of  the  somewhat  limited  capabilities  of  die  graphics  display  system  used  for  these  studies  and 
the  very  fine  arcs  required  to  represent  each  fine-range  resolution  cell,  a  conscious  decision  was  made 
to  display  range-azimuth  data  in  a  rectangular  format  to  prevent  any  loss  of  visual  information.  The 
drawback  to  this,  as  will  be  seen  in  a  following  chapter,  is  that  one  distorts  the  planar  geometry  and 
makes  comparisons  with  the  results  of  light  model  calculations  more  diffictilt  than  if  both  could  be 
displayed  in  range-azimuth  space.  Although  the  shape  of  the  data  display  in  the  rectangular  fonnat  is 
not  going  to  be  correct,  an  effort  was  made  to  make  the  width/depth  ratio  for  the  displayed  data  the 
same  as  the  azimuth/range  depth  ratio  in  the  real  data. 

Another  concern  about  how  to  best  visualize  the  data  and  calculation  results  for  each  test  is  how  to 
m^e  certain  that  one  test  can  be  immediately  compared  with  another.  In  other  words,  nothing  done 
during  the  compression  of  data  or  the  application  of  metrics  to  those  data  should  change  the  relative 
magnitudes  of  the  results.  One  should  be  able  to  look  at  a  color  display  of  power  returns  for  test 
MS5502  and  one  for  test  MS6701  and  be  confident  that  what  appears  as  a  given  color  on  one  test  is 
comparable  with  the  same  color  on  another  test.  To  this  end,  some  effort  was  made  to  determine  the 
bounds  of  data  and  metric  applications  to  all  of  the  data  before  selecting  a  color  bar.  The  absolute 
magnitudes  of  test  data  are  meaningless  because  the  data  are  uncalibrated. 
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An  example  of  the  data  is  shown  on  Figure  3,  which  contains  three  different  representations  o!  the 
measure  of  power  relumed  to  the  radar  from  the  terrain  and  targets  at  the  site  for  test  MS6603.  On  a 
linear  scale,  only  a  few  features  appear  as  bright  colors  representing  strong  returns.  As  will  be  seen 
later,  however,  these  are  not  all  man-made  targets.  This  type  of  display  is  a  kind  of  data  thresholding, 
and  is  very  much  a  function  of  the  strongest  signal  within  the  scene.  It  would  be  a  simple  matter  to 
defeat  a  sensor  system  based  simply  on  the  magnitude  of  power  returned  to  a  radar  through  the  use  of 
radar  comer  reflectors  that  arc  optimized  for  the  wavelength  in  question.  The  accompanying  three- 
dimensional  representation  of  the  linear  data  shows  just  how  complex  the  terrain  returns  can  be  for  this 
test  site.  It  is  not  at  all  apparent  that  there  are  four  man-made  targets  within  this  patch  of  natural 
terrain,  three  of  which  are  totally  in  the  open  and  one  is  camouflaged  with  brush.  More  will  be  said 
about  this  observation  later.  Terrain  elevation  data  contained  in  Appendix  C  reveal  that  the  low 
returns  at  the  longest  slant  ranges  are  due  to  shadowing. 

The  preferred  representation  of  data  in  the  range-azimuth  plane  is  logarithmic  such  as  that  shown  in 
part  (c)  of  Figure  3.  Here  it  is  apparent  that  there  are  several  areas  of  strong  returns  and  others  that 
are  either  in  shadow  or  that  absorb  most  of  the  millimeter  wave  energy.  The  black  ovals  represent  the 
approximate  locations  of  target  vehicles  within  this  scene.  Although  many  target  detection  algorithms 
probably  operate  in  the  linear  domain,  there  appears  to  be  much  more  informational  content  in  a  loga¬ 
rithmic  display,  particularly  when  one  is  trying  to  relate  terrain  features  to  measured  data  and  the 
results  of  calculations.  Note  in  the  two-dimensional  displays  that  a  histogram  of  data  values  has  been 
included.  These  will  be  significant  in  determining  the  quantitative  relationships  among  various  tests 
and  will  be  further  discussed  in  the  next  chapter. 
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MEASURE  OF  POWER 


MS6603_BAR0  .  +  LR**e 
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a.  Slant  range-azimuth  plot,  linear  scale 


Figure  3.  Measure  of  power  returned  from  test  MS6603  (Sheet  1  of  3) 
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b.  Three-dimensional  representation,  x  =  azimuth,  y  =  slant  range,  z  =  measure  of  power  (linear 
scale) 

Figure  3.  (Sheet  2  of  3) 
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3  Metrics  Approach  To  Quantifying 
Target-Like  Features 


Philosophy 

Aided  Target  Recognizer  (ATR)  algorithms  for  modem  sensor  systems  often  begin  the  process  of 
locating  and  identifying  targets  in  cluttered  backgrounds  by  eliminating  substantial  portions  of  the  data 
that  are  not  likely  to  contain  real  targets.  This  process,  referred  to  as  segmentation,  usually  involves 
the  application  of  a  gross  (in  a  spatial  sense)  signal-to-clutter  filter  to  coarse  range  profiles  and 
assumes  that  man-made  metallic  targets  produce  strong  returns  compared  with  the  natural  backgrounds. 
Only  after  areas  having  high  probabilities  of  containing  a  target(s)  have  been  identified  do  the  algo¬ 
rithms  then  priKAJss  data  ai  its  finest  spatial  resolution  to  complete  the  search  and  identification  proce¬ 
dure.  To  answer  questions  about  the  ability  of  natural  target  backgrounds  to  produce  target-like 
features,  one  must  not  throw  away  any  data.  Thus  the  approach  taken  in  these  studies  is  to  generate 
and  retain  the  finest  spatial  resolution  information  that  is  available  from  these  data. 

It  would  be  an  enormous  task  (as  well  as  an  infringement  on  the  proprietary  nature  of  many  sensor 
system  algorithms)  to  apply  all  of  the  target  feature  algorithms  to  these  data  that  are  currently  being 
used  or  being  considered  by  ATR  developers.  Therefore,  a  decision  was  made  to  apply  only  the  most 
commonly  used  and  simplest  features  as  a  first  step  in  tabulating  statistics  on  the  occurrence  of  target- 
like  signatures  in  the  MSFD  test  site  backgrounds. 

The  basic  approach,  then,  to  quantifying  the  occurrence  of  target-like  features  in  natural  back¬ 
grounds  is  to  first  select  several  simple,  widely  accepted  target  filters,  or  metrics,  for  further  study. 
These  metrics  are  then  applied  to  the  finest  resolution  range  profiles  for  a  representative  set  of  test  site 
data,  and  the  results  are  tabulated  and  stored  in  data  files  for  further  processing.  Color  graphics 
displays  of  the  results  are  generated  in  range-azimuth  space  for  qualitative  interpretation.  The  metric, 
or  filtered,  data  files  are  then  analyzed  using  statistical  routines,  to  determine  moments  of  the  histo¬ 
grams  and  other  relevant  information  such  as  the  relative  positions  of  target  data  within  the  filtered 
data  histograms.  This  latter  step  requires  placing  range  windows  over  the  known  target  locations  and 
eliminating  the  target-related  data  from  the  background  data  to  develop  background-only  histograms  or 
to  simply  quantify  the  target  signature  with  respect  to  the  background. 
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Target  Feature  Metrics 


Following  numerous  discussions  with  sensor  developers,  four  simple  target  feature  metrics  were 
selected  for  application  to  the  MSF^  fine-range  resolution  test  site  data  files.  These  are  listed  in 
Table  2  along  with  some  comments  about  their  general  applicability.  Even  though  the  amplitude  met¬ 
ric  does  not  require  any  additional  data  filtering,  it  proves  to  be  a  strong  indicator  of  target  pre.sence; 
as  long  as  ail  of  the  test  site  data  sets  are  properly  range  corrected,  amplitude  measurements  from 
different  sites  and  backgrounds  can  be  compared. 


Table  2 

Simple  Target  Feature  Metrics 

Metric  Oefirtition  Comments 


Amplitude  Measure  of  power  for  each  fine-range  Appropriate  tor  range  corrected  data.  M(»t 

resolution  cell.  useful  for  quick-look,  hot  spot  identification 

Signal-to-Clutter  (S/C)  Ratio  of  average  measure  of  power  in  a  Data  need  not  be  range  corrected  if  clutter 

number  of  fine-range  cells  forming  a  target  window  is  not  too  large.  Typical  applica- 
window  to  the  average  in  a  number  of  similar-  tions  are  region-of-mterest  identification  and 
size  clutter  windows.  One  buffer  window  is  target-size  S/C  for  detection, 
skipped  on  either  side  (in  range)  of  the  target 
window 

Variance  Variance  of  signal  (not  measure  of  power)  An  indicator  of  the  presence  of  one  or  more 

within  a  window  of  fine-range  resolution  cells.  strong  scatterers.  Data  need  not  be  range- 

corrected. 


Polarization  Ratio  Ratio  of  total  measure  of  power  in  the  left  Man-made  objects  should  exhibit  strong, 

receive,  right  transmit  channel  to  the  power  in  double  (even)  bounce  behavior.  Relatively 

the  right  receive,  right  transmit  channel.  smooth  terrain  should  have  strong,  single 

(odd)  bounce  characteristics.  Foliage 
should  be  neither.  Large  tree  trunks  and 
rock  outcroppings  are  expectec  to  cause 
problems. 


Tltere  was  some  concern  regarding  how  many  range-resolution  cells  should  be  included  in  each 
metric.  For  example,  the  signal-to-clutier  metric  consists  of  a  target  window,  a  buffer  window  on 
either  side  of  the  target  window,  and  a  number  of  clutter  windows  outside  of  the  buffers.  After  much 
experimenting  with  the  size  of  these  windows,  a  somewhat  meaningful  filter  size  was  chosen  that 
includes  a  2-m  target  window,  one  2-m  window  on  either  side,  and  four  2-m  clutter  windows  outside 
of  each  buffer.  Although  this  may  not  be  the  optimum  arrangement  from  a  sensor  developer’s  per¬ 
spective,  it  formed  a  reasonable  tradeoff  between  spatial  resolution  and  metric  amplitudes.  Larger  size 
windows  appeared  to  do  nothing  more  than  smear  the  results.  Other  filter  sizes  included  a  10-m  win¬ 
dow  for  variance  calculations  and  a  2-m  window  for  calculating  the  ratio  of  even  to  odd  polarization 
returns. 
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The  procedure  for  applying  each  target  Icaiure  metric  tcxccpi  for  amplitude)  was  to  sot  the  filter  at 
the  shortest  range  for  each  a/imuth,  calculate  the  value  ol  the  mcinc.  and  insert  dial  value  into  tlic 
fine-resolution  cell  at  the  middle  of  the  filler.  The  filter  was  then  moved  out  in  range  by  one  flnc- 
resolulion  cell  and  the  calculation  pcrlbmicd  again.  This  was  repeated  for  each  a/imuth  until  the  filter 
"bumped"  into  the  farthest  range  cell.  In  this  way,  some  data  is  lost  at  the  near-  and  far-range  edges 
of  the  range-azimuth  data  set. 


Representative  Results 

Range-azimuth  displays 

Figures  4-7  contain  graphics  display  results  of  the  application  of  each  target  feature  metric  to  the 
same  test  discussed  previously,  MS6603.  Each  figure  contains  the  slant  range-azimuth  representation 
in  both  linear  and  logarithmic  scales  along  with  a  crude  histogram.  Open  circles  added  to  the  logarith¬ 
mic  displays  approximately  coincide  with  the  locations  of  man-made  targets  at  the  test  site. 

Several  observations  are  in  order  at  this  time. 

a.  Of  the  four  different  representations  of  target  feature  metrics,  only  the  polarization  metric 
displays  a  log-normal  distribution  of  terrain  signatures.  Power  returns  from  natural  terrain  are 
often  assumed  to  be  distributed  log-normally,  but  the  power  returned  from  this  test  site  is  most 
assuredly  not  log-normal.  This  observation,  which  is  generally  true  for  the  other  test  sites,  is 
probably  due  to  a  combination  of  low-grazing  angles,  relatively  sparse  vegetation,  and  the 
presence  of  small  hills  and  depressions  that  result  in  numerous  "shadows"  within  each  test  area. 

b.  Returned  signal  strength  and  the  amplitudes  of  filtered  data  alone  are  not  enough  to  ensure  the 
detection  of  man-made  targets.  An  examination  of  all  of  the  MSFD  data  revels  that  if  one 
knows  where  these  targets  are  located,  the  data  will  often  show  a  strong  signature  at  about  that 
point  in  ^ace;  but  oUier  numerous  strong  features  appear  to  be  caused  by  either  the  terrain  or 
by  man-made  clutter  within  each  test  site.  As  for  man-made  clutter,  a  consensus  exists  among 
MSFD  participants  that  all  clutter  was  not  adequately  logged,  and  this  study  was  initiated  too 
long  after  the  demonstration  to  gather  any  reliable  information  by  revisiting  the  site. 

c.  The  polarization  ratio  concept  is  of  no  value  whatsoever  as  a  means  of  separating  man-made 
objects  from  natural  backgrounds.  After  some  reflection  on  this  question,  it  appears  that  even  if 
one  had  an  ideal  target-size  double-bounce  objea  within  the  radar’s  instantaneous  field  of  view, 
its  effect  on  the  returned  signal  would  still  be  minimized  because  of  the  fact  that  the  volume  in 
space  occupied  by  a  fine-range  cell  and  the  half  power  beam  width  of  the  radar  is  large  enough 
to  result  in  a  volume  averaged  response  that  is  essentially  depolarized.  For  this  test  case 
example,  the  half  power  beam  width  covers  an  arc  with  nominal  length  of  over  30  m.  The 
same  argument  could  be  made  to  account  for  why  so-called  target  amplitudes  are  not  greatly 
above  those  of  the  natural  background  in  many  test  cases. 
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Figure  4.  Range-azimuth  displays  of  the  ’measure  of  power"  metric,  MS6603 
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b.  Logarithmic  scale 

Figure  5.  Range-azimuth  displays  of  the  "signal-to-clutter"  metric,  MS6603 
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b.  Logarithmic  scale 

Figure  6.  Range-azimuth  displays  of  the  "variance"  metric.  MS6603 
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Figure  7.  Range-azimuth  displays  of  the  "polarization"  metric,  MS6603 
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(i.  While  all  of  the  targci  fcaiurc  metric  examples  shown  arc  Ibr  MSFD  data  that  are  derived  from 
sums  of  squares  of  measured  voltages  for  the  two  different  p<')larizaiion  combinations,  the  statis¬ 
tics  will  .show  that  neither  of  the  polarization  combinations  by  ihemvelves  consistently  per- 
fornied  better  than  the  other  when  looked  at  from  the  target  feature  metrics  perspective.  This  is 
a  reflection  of  observation  c. 


Range  profiles 

CotKider  now  another  visualization  of  these  target  feature  metrics,  this  time  as  a  trace  of  the  metrics 
along  a  given  azimuth  mark.  Several  of  these  range  profiles  are  shown  in  Figures  8-10.  It  is  obvious 
from  the  logarithmic  representation  of  power  for  test  MS6603  shown  in  Figure  4  that  a  man-made 
target  should  be  quite  apparent  along  azimuth  mark  18  (the  18th  vertical  strip  counting  from  the  left); 
those  profiles  are  shown  on  Figure  8.  The  apparent  target  position  is  indicated  by  the  dashed  line. 
Another  target  should  exist  along  azimuth  mark  21,  but  it  af^ars  to  be  obscured  by  strong  foreground 
returns.  Azimuth  mark  21  profiles  are  shown  on  Figure  9.  Finally,  there  are  numerous  significant 
returns  from  azimuth  mark  24  where  there  should  be  no  targets  of  interest.  Its  profiles  are  shown  on 
Figure  10. 

As  with  the  two-dimensional  maps,  several  observations  are  in  order  relevant  to  these  one-dimen¬ 
sional  traces. 

a.  With  respect  to  the  target  in  azimuth  mark  18,  its  measure  of  power  does  exceed  all  other 
returns  along  that  profile.  The  other  significant  return  at  about  range  cell  530  is  simply  a  spill¬ 
over  fiom  the  target  located  in  azimuth  mark  17.  At  this  point  it  is  impossible  to  speculate  as 
to  the  significance  of  the  shape  of  this  profile  because  not  enough  ground  truth  data  exist  to 
indicate  whether  the  reduced  returns  on  either  side  of  the  target  are  due  to  the  shadowing  effects 
from  adjacent  targets  or  dips  and  rises  in  the  local  elevation,  or  whether  this  is  Just  a  strong 
target  signature  superimposed  on  a  flat  terrain  return.  A  closer  examination  of  the  target  profile 
would  have  to  deal  with  higher  order  target  detection  algorithms,  anyway,  and  that  was  not  the 
objective  of  these  studies. 

b.  Regardless  of  what  influences  the  power  profile  data,  ^plication  of  the  signal-to-clutmr  metric 
to  a  strong  isolated  return  has  to  produce  the  characteristic  profile  shown  on  pan  (b)  of  the 
figure.  For  this  target,  under  these  test  conditions,  nothing  in  the  background  dominates  the 
signal-to-clutter  space  quite  like  the  real  target.  There  is  one  minor  anomaly  at  about  range  cell 
900  that  requires  more  ground  truth  than  that  available  to  explain. 

c.  By  its  very  definition,  a  calculadon  of  the  variance  of  a  signal  within  some  window  is  a 
measure  of  how  extreme  the  variations  in  the  signal  are.  One  is  therefore  not  surprised  to  see 
the  strong  signature  in  variance  space  because  of  this  rather  isolated  and  dominant  target  Large 
shadows  will  cause  the  severe  depression  of  the  variance  metric  seen  between  range  cells  900 
and  1100. 

d.  Even  for  the  well-defined  target  in  azimuth  mark  18,  the  polarization  metric  did  not  reveal  any 
anomalies,  as  expected. 
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Figures.  (Concluded) 
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Figure  10.  Range  profiles  for  target  metrics,  azimuth  mark  24  (Continued) 
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e.  Whereas  the  strong  target  signature  on  azimuth  mark  18  was  vcr\  isolated,  the  target  location 
on  azimuth  mark  21  docs  not  stand  out  in  any  of  the  target  feature  metric  profiles.  The  test 
summary  indicates  that  this  target  should  have  been  in  the  open  but  was  turned  at  an  angle  to 
the  radar.  In  other  words,  target  signatures  arc  not  always  the  strongest  signatures  on  any  given 
profile.  In  fact,  almost  everything  in  the  foreground  of  this  target  gave  stronger  returns.  It  will 
be  shown  later  that  the  strong  foreground  signals  can  be  attributed  to  terrain  geometry  and 
vegetation. 

/.  Finally,  consider  the  target  feature  metric  profiles  shown  on  Figure  10  for  azimuth  mark  24. 
Without  further  ground  truth  information,  it  appears  that  there  are  no  man-made  targets  along 
this  profile,  but,  at  least  in  the  signal-io-cluiter  domain,  there  are  strong  indications  of  target- 
like  features  between  range  cells  500  and  900.  The  very  elevated  power  measurements  within 
the  first  400  range  cells  are  uniformly  high  and  produce  no  target-like  features.  It  appeai:s  that 
anything  within  the  natural  terrain  that  could  cause  a  short  range  jump  in  power  return  could 
pn^uce  a  target-like  signal-to-clutter  response.  Within  the  terrain  at  Fort  Hunter  Liggett,  such 
behavior  could  be  the  result  of  the  isolated  trees  at  each  test  site. 


Statistics 

During  early  discussions  with  the  AATD  staff,  it  was  suggested  that  by  taking  a  hard  look  at  a  good 
set  of  existing  data  and  applying  the  metrics  analysis  approach,  a  quantitative  way  might  be  developed 
to  characterize  a  given  test  site  within  a  given  portion  of  the  spectrum.  After  a  great  deal  of  effort 
was  expended  on  this  task,  the  conclusion  to  date  is  that  a  simple,  clean  way  of  quantifying  a  test  site 
through  the  "eyes"  of  a  given  radar  system  has  not  been  found.  There  are,  however,  hints  that  such  an 
approach  might  still  woik.  Ideally,  what  would  be  required  is  well-documented  data  from  other  geo- 
gnqjhical  locations  that  were  collected  using  the  same  sensor  system.  Such  data  do  not  exist.  Never¬ 
theless,  fhe  following  observations  derived  from  statistical  analyses  may  prove  useful  for  future 
measurement  and  analysis  programs. 

Global  target  prominence.  Literature  on  imaging  ATR  systems'  has  defined  a  metric  concept 
which  quantifies  the  conspicuity  of  a  target  portion  of  an  image  relative  to  the  entire  image.  This  is 
referred  to  as  a  "global  target  prominence"  (GTP)  and  is  measured  by  determining  the  percentile  loca¬ 
tion  of  a  target  feature  value  relative  to  the  histogram  of  this  feature  for  the  entire  image  or  for  the 
poition  of  the  image  devoid  of  targets.  This  is  illustrated  in  Figure  11.  GTP  is  interpreted  as  the 
portion  of  the  image  less  conspicuous  than  the  target  for  a  specific  feature.  Values  range  between  0.0 
and  1.0.  A  value  of  1.0  indicates  that  the  target  has  a  higher  feature  value  than  anything  in  the  image; 
conversely,  a  value  of  0.0  indicates  that  the  target  has  a  lower  value  than  anything  in  the  image.  This 
concept  is  directly  adaptable  to  the  radar  features  used  in  this  study. 

Table  3  contains  average  GTP  values  for  each  of  the  radar  target  feature  metrics.  The  "all  data” 
column  represents  the  average  of  target  prominences  when  the  target  data  is  considered  as  pan  of  the 
histogram  for  each  test  The  "backgrounds  only"  column  uses  histograms  for  which  the  target  signa¬ 
tures  have  been  masked.  What  these  numbers  say  is  that  the  targets,  in  general,  are  among  the  domi¬ 
nant  features  within  all  of  the  test  sites.  However,  given  the  number  of  fine-range  cells 


'  J.  Beard,  L.  Clark  and  V.  Velton.  (1985).  "Characterization  of  ATR  performance  in  relation  to  image  measurements," 
unpublished  pspet,  AFWAIVAARF,  Wright-Patterson  AFB,  Ohio. 

26 


Chapter  3  Metrics  Approach  To  Quantifying  Target-Like  Features 


value  of  feature  X 


Figure  1 1 .  Graphical  representation  of  global  target  prominence  (GTP) 


Table  3 

Average  Global  Target  Prominence  (GTP)  of  Radar  Features 


(Baaed  on  $3  Targeta  In  24  Scenes;  Total  Measure  of  Poteer  Data) 

Measure  of  Target-Like  Featuies 

All  Dau 

Backgrounds  Only 

Pov»or 

0.994 

0.994 

Signal-to-Chitler 

0.990 

0.990 

Variance 

0.973 

0.975 

Polarization  Ratio 

0.925 

0.925 

contained  within  each  test  site  data  set,  there  are  still  a  large  number  of  other  features  whose  signa¬ 
tures  are  equally  as  strong. 

Cuimilative  histogram  inferences.  One  of  the  most  promising  approaches  to  quantitatively 
characterizing  test  sites  with  respect  to  a  given  radar  sensor  was  to  look  for  predictable  patterns  in  the 
histograms  of  measured  and  filtered  data.  The  following  figures  and  paragraphs  (x>ntain  some  observa¬ 
tions  regarding  the  cumulative  histograms  of  measure-of-power  data  for  all  of  the  test  sites  contained 
within  this  study. 

Figure  12  contains  all  of  the  MSFD  test  site  cumulative  histograms  presented  on  a  logarithmic 
scale.  There  is  no  obvious  clustering  of  these  curves  that  could  be  attributed  to  something  unique 
about  the  terrain  at  different  sites  within  Fort  Hunter  Liggett.  In  fact,  disregarding  the  two  anomalous 
responses  from  tests  MS6701  and  MS5803,  all  of  the  remaining  test  results  for  power  measurements 
are  quite  closely  clustered.  There  are,  however,  a  couple  of  observations  worth  noting. 
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Figure  12.  Cumulative  histograms  of  power  data  (all  tests) 

Figure  13  is  another  plot  of  cumulative  histograms  for  test  sites  whose  average  range  is  different 
from  all  of  the  others.  In  a  gross  sense,  increasing  range  to  the  test  area  appears  to  cause  a  shift  to  the 
right  in  the  cumulative  histogram.  The  net  result  of  such  a  shift  is  to  decrease  the  extremes  of  contrast 
within  the  test  area  that  could,  in  turn,  lead  to  less  apparent  target  signatures.  All  of  this  could  be 
exjdained  by  the  large  range  cell  volume  averaging  concept  discussed  earlier.  When  one  combirtes  a 
strong  target  signature  with  a  larger  number  of  worker  background  signatures,  the  resulting  average 
tends  toward  the  weak  background  levels. 
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Figure  13.  Cumulative  histograms  showing  range  effect 

The  final  observation  in  this  section  on  statistical  perspectives  is  that  the  data  collected  at 
Fort  Hunter  Liggett  do  appear  to  be  quite  repeatable.  Figures  14  and  15  show  the  histograms  of  tests 
whose  taiget  arrays  were  locked  at  about  the  same  position  relative  to  the  radar  but  whose  data  were 
collected  at  different  times  of  the  day  varying  from  predawn  to  the  middle  of  the  afternoon.  This 
clearly  shows,  also,  that  environmental  factors  such  as  air  temperature  and  relative  humidity  changes 
within  reasonable  limits  have  little  effect  on  radar  measurements  taken  during  MSFD.  There  is  some 
temptation  to  attribute  the  differences  between  the  histograms  shown  in  Figures  14  and  15  to  gross 
differences  in  weather  conditions,  as  the  data  represented  on  Figure  14  was  collected  during  heavy  fog 
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Figure  14.  Cumulative  histograms  showing  repeatability  (heavy  tog) 
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Figure  15.  Cumulative  histograms  showing  repeatability  (moderate  weather  conditions) 
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conditions  and  those  shown  on  Figure  15  to  a  more  moderate  set  ol  relative  humidit\  conditions. 
However,  the  magnitudes  of  the  data  do  not  suppon  tltat  contention.  Water  will  tend  to  absorb  and 
forward  scatter  the  radar  energy.  If  moisture  was  the  controlling  tacior  in  the  dilfercncc  between  these 
two  groups  of  tests,  then  the  data  in  Figure  14  should  be  shifted  to  the  left  ol  the  data  in  Figure  15. 
and  that  is  not  the  case.  The  range  effects  noted  above  seem  to  be  the  best  argument  lor  the 
differences. 
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4  Backscatter  Modeling  and 
Discussions  of  Terrain  Effects 


A  Simple  Backscatter  Prediction  Model 

The  second  approach  taken  in  these  studies  of  characterizing  the  interaction  of  the  MSFD  test  site 
backgrounds  with  the  radar  sensors  was  to  develop  and  exercise  a  simple  backscatter  prediction  model 
to  qualitatively  assess  the  impact  of  terrain  geometry  on  measured  data.  In  other  words.  "How  much 
of  the  radar  return  can  be  pr^icted  by  some  reasonable  knowledge  of  the  test  site  surface  geometry?” 
This  was  thought  to  be  a  particularly  relevant  question  for  the  MSFD  because  of  the  relatively  clean 
environment  of  the  tests  compared  with  what  might  exist  at  test  sites  in  other  temperate  parts  of  this 
country  or  in  western  Europe,  or  even  Central  America. 

The  backscatter  model  developed  consisted  of  the  following  elements  that  are  schematically  drawn 
on  Figure  16.  First  of  all,  available  25-m  terrain  elevation  data  were  converted  into  a  three-dimen¬ 
sional  finite  element  surface  of  triangular  facets  for  which  vectors  nc.mal  to  each  facet  are  easily 
computed.  Next,  a  point  light  source  was  positioned  at  the  location  of  the  radar  system  (either  the 
hi^  site  or  the  low  site),  and  the  angle,  0,  between  the  light  source  and  the  normal  vector  to  each 
terrain  facet  was  computed.  Under  the  assumption  of  Lambertian  scattering,  the  intensity  (or  power) 
of  tire  backscattered  light  was  calculated  to  be  proportional  to  the  cosine  of  this  angle.'  A  ray-tracing 
scheme  was  developed  to  ensure  that  terrain  facets  that  could  not  be  seen  by  the  radar  would  not  pro¬ 
duce  a  returned  signal.  Because  none  of  the  avaOable  Geographic  Information  System  packages  avail¬ 
able  at  WES  possessed  a  finite  range  point  light  source  simulation  capability,  this  software  package 
was  developed  by  the  first  author. 

The  point  light  source  simulation  of  a  radar  is  reasonable,  because  light  energy  and  radar  energy  are 
the  same  physical  phenomenon;  they  differ  only  in  frequency.  The  Lambertian  assumption  for  scatter¬ 
ing  is  not  unreasonable,  because  at  the  frequency  of  operation  of  this  radar  system,  the  terrain  surface 
is  certainly  rough.  Ultimately,  the  acid  test  for  this  model  is  its  comparison  with  real  data.  Results 
indicate  a  strong  correlation  between  the  two. 

Representative  results  for  this  simple  backscatter  prediction  model  are  shown  on  Figures  16-18. 

First  of  all,  Fig'.tres  16  and  17  show  the  terrain  backscatter  predictions  for  the  entire  areas  widiin 
which  MSFD  tests  were  conducted,  with  Figure  16  depicting  the  radar  located  at  the  high  site  (Site  8) 


‘  M.  Bom  and  E.  Wolf.  (2980).  Principles  of  optics.  6th  Edibon  (with  concctions).  Pergamon  Press,  Inc.,  Oxford.  UK. 
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Figure  16.  Schematic  of  the  backscatter  model 

and  Figure  17  showing  results  for  the  radar  located  at  the  low  site.  There  is  a  black  dot  on  each  fig¬ 
ure  at  the  location  of  the  radar,  and  40-ft  contours  are  superimposed  on  the  figures.  The  simple  model 
predicts  intensities  that  vary  between  0.0  and  1.0  in  magnitude.  These  results  were  then  converted  to  a 
40-db  range  (which  is  comparable  witti  the  range  of  measured  data)  and  displayed  with  the  same  color 
bar  used  for  the  displays  of  target  feature  metric  calculations  shown  in  the  previous  chapter.  Even  on 
the  gross  scale  of  these  figures,  it  is  readily  apparent  that  the  model  behaves  correctly  in  the  sense  that 
steep  slopes  visible  to  the  radar  cause  strong  backscatter  returns  and  that  tenain  that  caruiot  be  seen  by 
the  radar  is  placed  in  shadows  (white  background  on  these  figures  to  facilitate  viewing  of  dte  elevation 
contours). 

But  the  real  question  is  whether  or  not  the  model  can  qualitatively  simulate  the  measured  data  at 
each  test  site,  thereby  establishing  that  surface  geometry  is  a  major  contributor  to  the  test  results  at 
MSFD.  Returning  orKe  again  to  the  same  test  for  which  results  have  already  been  discussed  in 
previous  chiq}ters.  Figure  18  contains  a  comparison  between  measured  data  and  predicted  results  using 
the  simple  backscatter  model.  When  viewing  these  images,  keep  in  mind  that  the  measured  data  were 
displayed  on  a  rectangular  format,  while  the  outline  on  the  simulation  results  depicts  the  polar  nature 
of  the  radar  operation.  A  similar  set  of  figures  is  shown  in  Appendix  A  for  each  MSFD  test.  They  all 
demonstrate  a  strong  correlation  between  pnediaions  based  on  simple  surface  geometry  and  reflectance 
models  and  the  actual  data  measured  by  the  radar.  For  example,  note  that  on  Figure  18  strong  returns 
were  measured  by  the  radar  and  displayed  on  the  lower  right  hand  area  of  the  range-azimuth  plot  that 
could  rx>t  be  attributed  to  the  presetKe  of  any  man-made  objects.  The  simple  backscatter  model  pre¬ 
dicted  the  same  area  of  strong  returns  because  of  the  fact  that  the  terrain  surface  is  neatly  normal  to 
the  radar  at  that  location.  In  fact,  most  of  the  signiricant  features  of  the  test  site  radar  response  are 


34 


Chapter  4  Backscatter  Modeling  and  Discussions  of  Terrain  Effects 


Figure  17.  Backscatter  model  predictions  for  the  high  site  location 

qualitatively  simulated  by  this  simple  reflectance  model.  Hill  ^opes  facing  the  radar  are  illuminated, 
and  shadows  exist  where  depressions  in  the  terrain  have  been  measured. 


Vegetation  Overlays 

Also  shown  on  Figure  18  is  a  lO-ft  contour  map  with  a  vegetation  overlay.  This  was  produced  by 
enlarging  a  set  of  serial  photographs  of  the  MSFD  test  areas  taken  within  the  next  year  after  the 
demonstration  was  conducted  to  the  same  scale  as  the  measured  data  and  prediction  displays.  A  con¬ 
tour  map  was  generated  on  a  transparency  and  laid  over  the  aerial  photo.  As  can  be  seen  by  the  site 
photographs  contained  in  Appendix  A,  the  crowns  of  the  trees  and  other  significant  vegetation  arc 
easily  recognized  and  were  simply  colored  in  on  the  transparency.  Such  an  overlay  would  probably  be 
most  useful  for  those  test  sites  where  the  terrain  geometry  was  not  such  a  strong  contributor  to  the 
backscatter  signal  as  is  true  for  MS6603. 
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Figure  18.  Backscatter  model  predictions  for  the  low  site  location 
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5  Summary  and  Recommendations 


In  an  attempt  to  develop  a  scene  characterization  methodology  for  long-range  (2  to  7  km),  low- 
grazing  angle  (less  than  10  deg)  active  millimeter  wave  radar  systems,  a  twofold  approach  was  taken. 
First,  slant  range  versus  azimuth  maps  of  radar  backscatter  measurements  were  filtered  using  the  target 
feature  metric  technique  to  collect  statistics  on  the  occurrence  of  target-like  features  and  to  relate  those 
signatures  to  natural  terrain  conditions.  Second,  a  prediction  of  radar  returns  from  the  terrain  was 
modeled  by  a  point  light  source  at  the  radar  location  and  the  assumption  of  Lambertian  scattering  from 
the  terrain  facets.  Model  predictions  were  combined  with  overlays  of  vegetation  at  each  test  site  and 
compared  qualitatively  with  the  measured  data  to  assess  the  impact  of  terrain  conditions  on  backscatter 
response. 

Application  of  these  analysis  tools  to  the  Martin  Marietta  radar  data  collected  during  the  MSFD 
resulted  in  the  following  observations; 

a.  Man-made  targets  within  the  terrain  backgrounds  found  at  Fort  Hunter  Liggett  are  reasonably 
easy  to  locate  in  range-azimuth  "images"  when  their  positions  are  known,  but  might  not  be 
easily  found  by  automated  methods  applied  to  either  raw  data  or  filtered  data  because  of  the 
strong  target-like  returns  from  the  terrain  at  some  test  sites. 

b.  Not  enough  data  on  different  types  of  target  backgrounds  are  available  from  the  MSFD  trials  to 
develop  an  effective  statistical  tool  for  quantifying  background  conditions  for  differing  terrain 
conditions.  The  use  of  cumulative  histograms  for  achieving  this  objective  is  still  a  potential 
tool,  but  more  data  from  different  backgrounds  and  collected  with  the  same  radar  system  is 
needed  to  pursue  this  train  of  thought. 

c.  A  simple  point  light  source  model  with  assumed  Lambertian  scattering  qualitatively  correlates 
very  well  with  observed  data.  Vegetation  overlays  serve  to  enhance  that  correlation.  This 
implies  that  for  the  type  of  terrain  offered  by  the  Fort  Hunter  Liggett  test  sites,  even  crude 
surface  geometry  models  and  vegetation  overlays  will  explain  much  of  the  background 
response. 

d.  Many  questions  about  the  impact  of  terrain  on  radar  returns  under  the  conditions  stated  above 
remain  unanswered  in  part  because  of  insufficient  supporting  data.  Future  tests  in  which  radar 
data  will  be  closely  studied  must  include  more  comprehensive  ground  truth  measurements.  For 
example,  had  it  been  known  prior  to  the  execution  of  the  MSFD  that  a  detailed  look  at  the  radar 
data  was  going  to  be  taken,  a  complete  set  of  low-altitude,  oblique,  high-resolution  color  pho¬ 
tography  and  video  of  all  of  the  test  sites  could  have  been  collected  (from  at  least  two  of  the 
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cardinal  directions).  This  would  have  answered  any  questions  about  vegetation  conditions  or 
about  the  occurrence  of  man-made  clutter  and  such  natural  features  as  rock  outcroppings. 

Having  developed  the  data  management  and  analysis  tools  for  this  study,  it  would  be  relatively  easy 
to  test  other  higher  order  target  feature  metrics.  In  fact,  a  database  now  exists  for  this  one  set  of  ter¬ 
rain  conditions  to  take  a  closer  look  at  existing  target  detection  algorithms  or  to  develop  and  test  new 
ones.  One  very  straightforward  task  that  could  be  undertaken  and  that  would  require  nothing  more 
than  to  revisit  these  existing  data  would  be  to  examine  the  impact  of  averaging  multiple  azimuth 
sweeps  or  comparing  results  from  different  elevation  settings  as  suggested  in  Chapter  2. 

In  the  area  of  modeling,  it  appears  from  the  simple  exercise  conducted  within  this  study  that  modi¬ 
fications  of  the  point  light  source  model  could  yield  a  valuable  tool  for  predicting  terrain  signatures. 
Examples  of  such  modifications  would  include  more  realistic  scattering  assumptions  from  grass- 
covered  terrain  and  a  coupling  of  tree  and  bush  spatial  information  with  a  separate  module  for  estimat¬ 
ing  an  average  backscatter  from  each  species. 
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Appendix  A 

Scene  Metric  Statistics 


The  following  pages  contain  a  summary  of  common  statistics  for  the  log-scale  histograms  associated 
with  each  of  the  24  MSFD  radar  scenes.  There  is  a  table  for  data  that  represent  the  total  power 
returned  to  the  radar  receiver  (BAND=total  power)  as  well  as  a  table  for  the  power  measured  on  the 
right  receive-right  transmit  channel  (BAND=polarization  RR)  and  for  power  measured  on  the  left 
rcceive-right  transmit  channel  (BAND=polarization  LR).  Each  of  these  tables  has,  in  turn,  separate 
sheets  for  each  radar  target  feature  metric  (METRIC=measure  of  power,  METRIC=signal  to  clutter, 
METRlC=variance).  The  final  table  summarizes  results  for  the  polarization  ratio  metric 
(BAND=LR/RR,  METRIC=polarization  ratio). 
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Target  Metric  Statistics 


In  an  attempt  to  make  this  document  as  useful  as  possible,  statistics  on  radar  feature  metrics  as  they 
relate  to  the  actual  targets  were  compiled  and  are  included  in  this  appendix.  The  procedure  for  com¬ 
piling  these  statistics  was  the  following.  First  of  all,  a  mask  was  created  for  each  scene  that  served 
two  purposes.  One  was  to  remove  the  influence  of  a  man-made  target  on  the  radar  return  by  effec¬ 
tively  removing  from  the  scene  a  portion  of  one  or  more  azimuth  slices  that  was  centered  on  the  likely 
location  of  each  target.  The  other  purpose  was  to  identify  the  strongest  return  from  that  segment  as 
the  target  return,  or  metric  value.  The  target  metric  value  was  then  compared  with  the  target- free 
background  to  produce  the  global  target  prominence  for  that  target  that  is  defined  in  Chapter  3  of  the 
main  text. 

As  with  the  previous  appendix,  data  are  presented  in  terms  of  both  the  type  of  received  signal  con¬ 
sidered  (total  power,  right  reccive-iight  transmit  (RR)  polarization  power,  left  rcceive-right  transmit 
(LR)  polarization  power)  and  the  metric  being  considered  (measure  of  power,  signal -to-clutter,  vari¬ 
ance,  and  polarization  ratio). 
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Appendix  C 

Relevant  MSFD  Test  Site  Data  and 
Analysis  Results 


Contained  within  this  appendix  is  all  of  the  infonnation  used  to  study  the  radar  rentms  from  each  of 
the  24  MSFD  test  sites  for  which  good  data  were  available.  Arranged  in  numerical  order  by  test  num¬ 
ber,  each  packet  of  infonnation  contains  a  test  summary  sheet,  Charge  Coupled  Device  B&W  oblique 
photographs  (if  available),  a  B&W  down-looking  aerial  photo  and  a  10-ft  contour  map,  both  two- 
dimensional  and  three  dimensional  {3-D)  linear  scale  representations  of  total  measured  power,  slant 
range-azimuth  representations  of  measured  data  and  filtered  data  from  the  target  feature  metric  exer¬ 
cises,  range-azimuth  representations  of  backscatter  modeling  results,  and  another  site  contour  map  with 
vegetation  overlays.  The  af^roximate  locations  of  targets  for  each  test  are  shown  as  circles  on  the 
backscatter  prediction  results. 
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Figure  C2.  Summary  data  for  MS5503  (Sheet  1  of  7) 
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(Sheet  3  of  6) 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


Figure  C3.  (Sheet  4  of  6) 
C18 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Resulte 


Figure  C3  (Sheet  6  of  6) 


C20 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


15  Mar  1988  0825  Test  ^?^MS5601 


Figure  C4.  Summary  data  for  MS5601  (Sheet  1  of  7) 


Figure  C4.  (Sheet  3  of  7) 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


C23 


Figure  C4.  (Sheet  4  of  7) 
C24 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


Figure  C4.  (Sheet  5  of  7) 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


C25 


Figure  C4.  (Sheet  6  of  7) 
C26 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


Appendix  C  Retovant  MSFD  Test  Site  Data  and  Analysts  Results 


15  Mar  1988  0905  Test  #MS5602 


q: 

<: 

s: 

u 

or. 


z 

IxJ 

CE 

D 


A  O 

03 


SSi 

fi.  0. 

o  o 


11 


o  a 

aCBsSas 
E 


Si 


iliisii 


□ 


c 


> 

LJ 


2£i*SS0S 

sssssss 


_1 

u 


X 

I— 

X 


N 

C 


$228898 

2*'98*'8$ 

M  A  M  i*  M  •« 


(/O 

LiJ 

Y- 

<Z 

X 

»— I 

a 

q: 

a 

a 

CJ 


u 

Q. 

>- 


iiiissa 


u 

o 

O' 

< 


s 

tt 


C28 


Appencfx  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


Figure  C5.  Summary  data  for  MS5602  (Sheet  1  of  7) 
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Measured  data,  linear  scale 


3-D  plot  of  measured  data,  linear  scale 
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Figure  C6  Summary  data  lor  MS5603  (Sheet  1  of  7) 
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Figure  C7.  Summary  data  for  MS5604  (Sheet  1  of  7) 
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Measured  data,  linear  scale 
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3'D  plot  of  measured  data,  linear  scale 
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Figure  C8.  Summary  data  for  MS5702  (Sheet  1  of  6) 
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Figure  C9.  Summary  data  for  MS5703  (Sheet  1  of  6) 
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Figure  CIO.  Summary  data  for  MS5704  (Sheet  1  of  6) 
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Measured  data,  linear  scale 


3-D  plot  of  measured  data,  linear  scale 
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Figure  C1 1 .  Summary  data  for  MS5801  (Sheet  1  of  6) 
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Measured  data,  log  scale 


Backscatter  predictions,  log  scale 
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Measured  data,  linear  scale 


3-0  plot  of  measured  data,  linear  scale 
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Figure  C13.  Summary  data  for  MS5803  (Sheet  1  of  6) 
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Figure  C14.  Summary  data  for  MS5804  (Sheet  1  of  6) 
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Measured  data,  linear  scale 


3-D  plot  of  measured  data,  linear  scale 
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Measured  data,  linear  scale 


Figure  C15. 
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3-D  plot  of  measured  data,  linear  scale 
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Measured  data,  linear  scale 
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3-D  plot  of  measured  data,  linear  scale 
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Measured  data,  linear  scale 
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Measured  data,  linear  scale 
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Measured  data,  linear  scale 
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Figure  C18  (Sheet  5  of  7) 
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Overhead  photo 
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Figure  C20.  Summary  data  for  MS6301  Short  Range  (Sheet  1  of  7) 


Figure  C20  (Sheet  2  of  7) 
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Measured  data,  linear  scale 


3-D  plot  of  measured  data,  linear  scale 
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Figure  C21 .  Summary  data  for  MS6301  Long  Range  (Sheet  1  of  7) 


Figure  C21.  (Sheet  2  of  7) 
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Overhead  photo 


Terrain  contours  (10-ft  interval) 


Figure  C21 .  (Sheet  3  of  7) 

Appendix  C  Relevant  MSFO  Test  Site  Data  and  Analysis  Results 

C135 

3-D  plot  of  measured  data,  linear  scale 

Figure  C21  (Sheet  4  of  7) 


Cl  36 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


Figure  C21  (Sheet  5  of  7) 

Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


Figure  C21  (Sheet  6  of  7) 
Cl  38 


Appendix  C  Relevani  MSFD  Test  Site  Data  and  Analysis  Results 


Figure  C21  (Sheet  7  of  7) 

Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


C139 


24  Mar  1988  0557  Test  ?5^MS6403 


Cl  40 


CO 


Wl 

ii 

IP 

fMSi 

ii 

H 

iWi 

£ 

m 

^21' 

ii 

to 

Ck; 

< 

2: 

UJ 


Z 

UJ 

a 


z 

a 


< 

> 

u 


X 

h- 

3 

Z 
»— < 
M 
C 


(/) 
UJ 
J- 
9  < 

■“  z 

h—l 

a 
a: 
□ 

o 


ui 
5>  z 


< 


s  ^ 

*  D. 
>- 


U 

a 

a: 

< 


8 


iiil 


0.  ».  A.  ft.  ^ 


III 


m 


SS;X$ 
99  ;;8 
S8SSS 


i&ssa 

5S!«2* 

SSSSI 


Hill 


isigg 


Appendix  C  Relevant  MSFD  Test  Site  Data  and  Analysis  Results 


Figure  C22.  Summary  data  for  MS6403  (Sheet  1  of  7) 


Figure  C22.  (Sheet  2  of  7) 
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Figure  C22.  (Sheet  7  of  7) 
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Figure  C23.  Summary  data  for  MS6603  (Sheet  1  of  7) 


Figure  C23.  (Sheet  2  of  7) 
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Figure  C24.  Summary  data  for  MS6701  (Sheet  1  of  7) 
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